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A  STEADY  STATE  LONGITUDINAL  MANPOWER  PLANNING  MODEL 


WITH  SEVERAL  CLASSES  BY  MANPOWER 
by 

Richard  C.  Grinold 

ABSTRACT 

In  the  Navy's  officer  system  there  are  several  classes  of 
manpower  (e.g.  pilots  and  submariners)  that  perform  specialised 
jobs  and  also  can  perform  several  types  of  non-special i zed  jobs 
(e.g.  military  training,  personnel  management,  etc.).  This  report 
is  concerned  with  the  general  prob.1  m  of  allocating  the  different 
types  of  jobs  among  the  several  classes  of  manpower. 

The  report  describes  a  model  that  constructs  a  personnel 
inventory  by  rank  for  each  of  several  manpower  classes  (pilots, 
etc.)  and  then  allocates  those  people  to  the  specialized  and 
common  jobs  that  they  are  allowed  to  do. 

The  idea  of  the  model  is  to  allow  a  policy  maker  to  quickly 
reconcile  billet  requirements  with  the  reality  of  available  acces¬ 
sions,  job  sharing  targets  between  classes,  and  continuation  rates 
of  the  different  manpower  classes. 

Different  allocations  can  be  produced  either  by  assuming 
the  values  of  a  few  key  variables,  or  they  can  be  generated  using 
an  optimization  scheme  that  sets  "allowable"  percentage  errors. 

Pour  optimization  variants  based  on  this  idea  are  described.  The 
report  contains  some  typica1  data,  the  results  of  calculations, 
and  a  description  of  computer  programs  used  to  solve  the  problem. 
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A  STEADY  STATE  LONGITUDINAL  MANPOWER  PLANNING  MODEL 


WITH  SEVERAL  CLASSES  OF  MANPOWER 


by 

Richard  C.  Grinold 


1 .  Introduction 

This  paper  outlines  the  construction  of  a  manpower  plan¬ 
ning  model  for  a  system  in  which  each  of  several  classes  of  man¬ 
power  are  assigned  to  several  categories  of  jobs.  The  jobs  are 
either  specialized  in  that  only  one  type  manpower  can  perform 
the  job  or  general  in  that  more  than  one  manpower  type  is  quali¬ 
fied  to  do  the  job. 

Both  jobs  and  manpower  are  broken  down  by  experience 
level.  A  job  with  a  certain  experience  level  must  be  filled  with 
a  person  that  has  the  same  experience  level. 

The  model  assumes  a  steady  state.  In  each  year  the  job 
requirements  are  the  same.  In  each  year  the  same  number  of  people 
enter  the  system  (bottom  level  entry  only)  ,  the  same  number  leave 
the  system  from  each  experience  level,  and  the  same  number  move 
up  from  each  one  level  so  that  the  inventory  of  people  at  each 
exerpience  level  and  in  each  manpower  class  remains  the  same. 

There  are  several  items  represented  in  the  model: 
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(i)  The  requirements  for  jobs  by  experience  levels. 

(ii)  The  requirements  for  people  by  experience  levels. 

(iii)  The  job  sharing  targets  for  common  jobs;  the  fraction 
of  common  jobs  we  would  like  to  see  allocated  to  each 
manpower  class . 

(iv)  The  people  sharing  targets  for  manpower  classes;  the 
fraction  of  each  manpower  class  that  we  would  like  to 
see  allocated  to  each  type  of  job. 

(v)  The  rate  of  accessions  in  each  manpower  class. 

(vi)  The  retention  of  manpower  in  each  class;  i  .e .  the 

relationship  between  accessions  in  each  class  and  the 
inventory  of  people  by  experience  level  in  each  class. 

The  model  is  designed  to  study  the  interaction  of  these 
factors.  It  is  easy  to  see  from  the  list  above  that  items  (i) 
through  (iv)  are  interrelated  and  unless  care  is  taken  in  specify¬ 
ing  these  goals  they  will  be  inconsistent.  Also  items  (ii) ,  (v) 
and  (vi)  are  strongly  related  in  that  (v)  and  (vi)  determine  the 
inventory  of  people  in  each  class  at  each  experience  level.  Only 
in  the  most  fortunate  circumstances  will  this  be  consistent  with 
the  personnel  requirements  (ii). 

Items  (i)-(iv)  are  targets.  They  may  be  unrealistic  and 
inconsistent  but  they  can  be  set  by  the  manpower  planner.  Items 
(v)  and  (vi)  ,  however,  depend  to  some  extent  on  the  behavior  of 
personnel  in  the  system  (vi)  and  on  the  system's  ability  to 
attract  qualified  people  to  each  class  (v) .  The  reader 


should  keep  in  mind  this  caveat  and  be  prepared  for  shifts  in 
personnel  behavior  in  response  to  changes  in  other  system  parameters 
or  in  external  factors. 

The  model  is  a  laboratory  for  testing  the  relationship 
between  these  factors.  We  can  try  one  set  of  policies  and  look 
at  the  discrepancy  between  job  requirements  and  actual  require¬ 
ments  as  well  as  the  discrepancy  between  personnel  requirement 
and  the  actual  inventory  in  each  class. 

The  model  may  be  used  in  several  ways.  The  easiest  is  for 
the  planner  to  stipulate  a  policy  and  then  examine  its  effects. 

A  more  difficult  procedure  is  for  the  planner  to  stipulate  a 
range  of  policies  and  then  use  some  type  of  optimization  scheme 
to  select  a  policy  within  the  range.  This  paper  contains  four 
examples  of  this  type  of  optimization.  Each  is  based  on  the  notion 
of  a  penalty  function  the  measures  the  discrepancy  between  desired 
job  and  personnel  requirements  and  what  is  actually  provided. 

We  take  two  types  of  penalty  functions:  piecewise  linear  and 
quadratic,  and  we  examine  two  broad  sets  of  policies:  allocating 
people  to  jobs  using  the  job  sharing  rules  (  ( i i i )  above),  or 
allocating  people  to  jobs  using  the  people  sharing  rates  ((iv) 
above) .  This  gives  us  four  combinations  and  thus  four  distinct 
optimization  models.  The  piecewise  linear  penalty  functions  lead 
to  linear  programming  models  and  the  quadratic  penalty  functions 
to  the  minimization  of  a  quadratic  form  subject  to  linear  equality 
constraints . 
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The  model  is  motivated  by  a  study  of  the  U.  S.  Navy's 
officer  corps.  We  shall  carry  an  example  using  that  system 
throughout  the  text  in  order  to  illustrate  each  idea. 

The  model  is  based  on  a  longitudinal  manpower  flow  model. 

This  type  of  system  is  described  in  depth  in  Grinold  and  Marshall 
[1],  Chapter  3.  The  use  of  piecewise  linear  penalty  functions 
to  measure  the  discrepancy  between  actual  performance  and  stated 
objectives  is  commonly  called  "goal  programming."  This  idea 
has  been  extensively  developed  by  Charnes  and  Cooper;  a  good 
review  can  be  found  in  [3].  The  use  of  quadratic  penalty  functions 
is  quite  common  in  the  optimal  control  literature;  its  use  in 
more  behavioral  settings  was  pioneered  by  Holt,  Modigliani,  Muth 
and  Simon  in  [2]. 

The  paper  consists  of  several  short  sections  each  dedicated 
to  a  specific  point  and  most  illustrated  by  the  example  carried 
thoughout  the  text.  The  model's  structure  is  given  in  Sections  1-14. 
Section  15  is  a  review  that  gathers  all  the  definitions  presented 
to  that  point.  Sections  16  to  20  describe  alternative  ways  to 
choose  an  allocation,  the  idea  behind  our  use  of  penalty  functions, 
and  the  specific  construction  of  the  piecewise  linear  and  quad¬ 
ratic  penalty  functions. 

Appendix  A  describes  the  two  linear  programming  models 
that  arise  from  the  use  of  piecewise  linear  penalties  and 
Appendix  B  the  two  models  that  stem  from  the  use  of  quadratic 
penalties.  Appendix  C  contains  some  sample  solutions. 
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The  organization  reflects  the  relative  importance  of 
the  topics.  The  structure  of  the  model  is  the  most  important 
and  it  is  stressed  and  reviewed  in  Sections  2-14.  The  use  of 
optimization  to  select  an  allocation  is,  in  the  main,  merely  a 
device  to  circumvent  the  difficulty  of  having  a  wide  range  of 
policy  choice.  Once  the  optimization  rule  is  set  the  planner 
has  a  direct  route  from  policy  to  result.  It  is  the  variation 
in  results  that  comes  from  changes  in  policy  that  will  be  of 
most  interest  to  the  planner.  The  optimization  is  a  device  to 
help  make  that  connection. 

One  final  point  should  be  made  before  describing  the 
model's  structure.  This  is  an  aggregate  planning  model.  It  is 
intended  to  test  policies  that  will,  in  turn,  provide  a  foundation 
for  the  day-to-day  operating  of  the  system.  This  model  certainly 
will  not  show  us  how  to  operate  the  system.  No  model  can  answer 
all  questions  simultaneously  and  ours  is  no  exception. 


2 .  Manpower 


There  are  K  different  classes  of  manpower  indexed 
by  k  =  1,2,...,K.  The  classification  scheme  is,  of  course, 
directly  related  to  the  objectives  of  the  model  builder.  In 
general,  the  classification  should  be  fine  enough  to  capture 
the  important  substitution  possibilities  and  economical  in  avoid 
ing  the  listing  of  all  possibilities  in  a  futile  attempt  to  repl 
cate  reality.  In  our  example  there  are  five  manpower  types. 


Manpower 

Type 

Officer 

Designator 

Description 

1 

llOx 

Women 

2 

lllx 

(and  116x) 

Surface  warfare 

3 

112x 

(and  117x) 

Submarine  warfare 

4 

1 31x 
and 

1 39x 

Pilots 

5 

1 32x 
and 

Naval  flight  officers 

137x 

Table  2-1:  The  manpower  classes. 
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J.  Stages 

Rach  of  fiver's  career  is  broken  down  into  I -stages  indexed 
by  i  *  1,2,...,  1.  The  stage*  oan  be  defined  in  many  ways.  The 
simplest  is  by  period  of  service,*  if  we  track  officers  for  26 
years  then  the  index  i  would  go  from  one  to  26.  In  our  example 
the  stages  roughly  coincide  with  the  time  period  in  which  officers 
hold  a  certain  rank.  In  general  for  i  «  1,2,...,  I,  stage  i 
will  run  f tom  time  of  service  s._j  to  s^.  The  maximum  length 
of  service  is  Sj  -  S  and  s^  -  0. 

A  person  in  manpower  class  k  and  in  stage  i  of  their 
career  will  be  cal  lev!  a  type  (i,k)  I'otson. 


St  aue 

l.OS 

IVscrtpt  ion 

l 

0-2 

Kits ign 

2 

2-4 

l.ieutenant --J . d. 

\ 

4-u 

Lieut  enant 

4 

»»-  14 

l.ieutenant  Commander 

r\ 

M-l^ 

Commander 

6 

1^-26 

Capt  a  in 

Table  1- 1 .  Stages  or  experience  levels. 


4 .  Manpower  Flows 


We  assume  a  steady  state  model.  In  each  year  the  same 
number  of  people  will  enter  each  manpower  class,  and  the  inven¬ 
tory  of  people  in  each  year  of  service  and  in  each  manpower  class 
will  remain  constant. 

For  class  k  =  1,2,...,K  let  y^  be  the  number  of  people 
entering  manpower  class  k,  and  let  a(k,s)  be  the  fraction  of 
these  accessions  that  remain  in  the  system  for  s  years.  The 
index  s  runs  from  0  to  S  years  and  the  function  ct(k,s)  is 
decreasing  in  s.  The  function  a(k,s)  is  frequently  called  the 
survivor  curve . 

Stage  i  of  a  person's  career  runs  from  length  of  service 


s^_1  to  s ^  (sQ  =  0) .  The  number  of  people  in  stage  k  is 


„ikyk  where  „.k 


is  defined  by 


w 


ik 


/  a(k,s)ds  . 
si-l 


In  our  example,  let's  take  a  hypothetical  survivor  curve 
for  pilots  (class  4) .  We  follow  these  officers  for  26  years. 
Figure  4-1  shows  the  survivor  curve.  The  shaded  areas  are 
integrals  under  this  curve  for  the  duration  of  each  stage. 

The  calculation  of  w^  can  be  approximated  using  a 
discrete  form  of  the  survivor  curve. 

The  coefficients  w^  can  be  interpreted  as  the  amount 
of  time  a  person  in  manpower  class  k  is  expected  to  spend  in 
stage  i.  A  crude  way  to  view  this  is  to  say  that  w^  is  the 
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product  of  the  probability  of  reaching  stage  i,  and  the  length 
of  stage  i.  The  probability  will  be  ”  si-l*  ' 


i 

Wik 

Si  *  Si-1 

Prob .  of 
reaching 
stage  i 

s 

1 

1.97 

2 

0.985 

m 

2 

1.85 

2 

0.925 

H 

3 

2.50 

5 

0.5 

G 

4 

1.40 

5 

0.28 

E 

5 

1.22 

5 

0.25 

S 

6 

0.54 

7 

0.77 

Table  4-1.  The  expected  waiting  time. 


There  are  J  different  typ 
J.  In  our  example  we  c 


Job  No. 

Billet 

Code 

1 

1000 

2 

1050 

3 

1110 

and 

1160 

4 

1120 

and 

1170 

5 

1310 

and 

1390 

6 

1320 

and 

1370 

7 

1300 

Table  5-1 


6.  Billet  Requirements 

The  I  by  J  matrix  B  contains  the  billet  requirements 
Since  the  career  stage  is  a  measure  of  experience  (e.g.  length  of 
service  or  rank),  the  requirements  are  broken  down  by  stage. 

b^j  the  number  of  stage  i  people  need  for  jobs  type  j 

The  billet  requirements  matrix  for  our  example  is  shown  below. 

The  billets  will  be  identified  by  the  index  pair  (i,j). 
Thus  we  shall  speak  of  type  (i,j)  billets  and  type  (i,k)  people. 


JOBS 


1 

2 

3 

4 

5 

6 

7 

s 

1 

409 

0 

2095 

688 

961 

521 

0 

H 

2 

1008 

0 

1883 

734 

1998 

1212 

1 

H 

3 

1806 

378 

2080 

844 

3572 

1399 

6  36 

G 

4 

1495 

599 

1464 

810 

1780 

547 

726 

E 

5 

1031 

470 

926 

501 

678 

55 

902 

6 

592 

490 

321 

149 

0 

0 

387 

Table  6.  The  billet  requirements  matrix  B. 
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7.  Billet  Sharing  Array 

For  each  i  the  K  by  J  matrix  F1  gives  the  share 
(fraction)  of  the  type  (i,j)  billets  that  should  be  performed 
by  manpower  type  k.  The  elements  of  F1  are  denoted  f£j» 
they  are  nonnegative,  and  the  rows  sum  to  one,  fj^  *  1. 

If  f£j  >  0 ,  we  interpret  the  fraction  as  a  goal;  we  would 
like  that  fraction  of  the  (i,j)  billets  filled  from  manpower 
class  k.  However,  we  interpret  fjV  =  0  to  mean  that  manpower 
class  k  is  not  qualified  to  fill  billet  (i,j). 

There  is  also  possible  confusion  between  the  idea  of  job 
sharing  and  a  people  sharing  concept.  A  manpower  class  may  per¬ 
form  two  types  of  jobs:  the  jobs  for  which  it  is  uniquely 
qualified  and  other  jobs  that  are  shared  among  several  manpower 
classes.  Another  way  to  look  at  the  allocation  of  the  common 
jobs  is  to  stipulate  a  fraction  of  the  inventory  of  (i,k) 
people  that  should  be  assigned  to  job  j.  That  idea  will  be 
considered  in  Section  12. 

To  save  space,  we  shall  only  give  one  of  the  matrices 
F1  (for  i  =  3)  used  in  our  example. 


JOBS 


1  2  3  4  5  6  7 


M 

C 

a 

0.7 

— 

A 

N 

L 

ms 

0.15 

.5 

1 

P 

■ 

0 

A 

0.06 

.2 

1 

W 

H 

E 

S 

0.05 

.18 

1 

.  57 

R 

S 

B 

0.04 

.12 

1 

.43 

Table  7-1.  The  job  sharing  matrix  Fl  for  i  =  3. 
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jdwso*oz>3: 


8.  Target  Allocations 

The  number  of  (i,k)  people  desired  in  job  (i,j)  is 

given  by 


t* 


kj  kj  ij 


Recall  that  b..  is  the  (i,j)  billet  requirement  and  .  is 
13  KJ 

the  fraction  of  those  billets  to  be  filled  by  manpower  class  k. 

The  target  allocations  used  in  our  example  (for  i  =»  3) 
are  shown  in  Table  8-1. 


JOBS 


1 

2 

3 

4 

5 

6 

c 

D 

1264.2 

L 

9 

270.9 

189 

2080 

A 

B 

108.36 

75.6 

844 

S 

B 

90.3 

68.04 

3572 

362.52 

S 

5 

72.24 

45.36 

1399 

273.48 

Table  8-1.  The  target  allocation  t? 
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9.  Actual  Allocation 

Let  the  variable  give  the  actual  number  of  (i,k) 

people  assigned  to  billet  (i,j).  If  it  was  possible  we  would 
like  a£j  -  t^...  That  would  ensure  that  all  targets  are  met. 
That,  however,  is  usually  impossible.  Indeed,  reconciling 
the  a£j  and  tj^  is  the  purpose  of  this  model. 

The  variables  ajV  do  not  mean  that  a^  individuals 
from  class  k  are  locked  into  type  (i,j)  billets  when  they 
are  in  stage  i.  It  does  mean  that  at  any  time  a£_.  people 
from  class  k  are  filling  type  (i,j)  billets.  Any  particular 
(i,k)  person,  may  spend  stage  i  in  several  billets. 

If  f£j  *  0,  then  an  (i,k)  person  is  not  qualified  to 
fill  an  (i,j)  billet.  So  if  f£j  =  0,  we  shall  require  that 
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10 .  Conservation 


Recall  that  w^y^  is  the  number  of  type  (i,k)  people. 
The  allocation  must  satisfy  the  following  conservation 
relation 


wikyk 


for  all  (i ,k) 


This  says  that  all  manpower  fills  some  job. 
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IX.  Excess  and  Deficit 

One  of  our  main  objectives  in  constructing  this  model 
is  to  compare  the  allocation  and  target.  We  could  look  at 
all  the  discrepancies  a^k  -  t^.  That  would  be  very  general, 
but  it  may  not  be  comprehensible.  We  have  chosen  to  concen¬ 
trate  on  more  aggregate  measures  of  excess  and  deficit:  the 
actual  number  of  (i,k)  people  and  the  actual  number  of 
(i,j)  billets  filled  as  compared  with  the  targets.  Our 
notation  for  these  measures  is 


K 

I 


k=l 


the  number  of  poeple  assigned 
to  billet  (i,j) 


the  discrepancy  between  billet 
(i,j)  assignment  and  requirements 


z 


ik 


pik 


f  aj 

j  =  l 


J 

l 

j=l  R:1 


the  number  of  (i,k)  people 


the  total  requirement  for  type 
(i,k)  people 


z 


ik 


pik 


the  discrepancy  between  the  inventory 
of  (i,k)  people  and  requirements. 
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12.  People  Sharing 

Recall  that  we  defined  our  targets  from  the  viewpoint 
of  billet  sharing.  Given  any  allocation  we  can  calculate  the 
people  sharing  fractions.  Let  g1  be  the  fraction  of  type 

^  j 

(i,k)  people  that  are  assigned  to  job  j,  thus  billet  (i,j). 


The 


Jkj 


are  given  implicitly  by  the  equation 


akj  gkjZik 


A',  mentioned  above,  an  alternative  model  could  be  constructed 
using  the  people  sharing  fractions  g^  as  input  targets 
and  the  billet  sharing  fractions  as  outputs.  In  that 


and  the 


case  we  would  start  with  people  requirement  pik, 

1  and  define  target  allocations  by 


people  sharing  rules  gkj 


fckj  gkjpik" 

Then  we  can  calculate  b. .  from 

il 

bij  "  j-  gkjPik' 

In  any  case, for  consistency,  the  following  relations  should 

hold 
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(i)  9kjpik  ■  £kjbij 


(iil  k  9kj  *  1  p 


ik 


ib  tk^ 


uli)  X  fkj  ■  1  “>  bij  ■  X  ^ik 


To  assure  consistency  we  should  either  start  with  f 
and  b  and  then  calculate  g  and  p,  or  start  with  g  and 
p  and  calculate  f  and  b.  If  we  have  a  desire  to  start  with 
g  and  b,  then  one  suggestion  is  to  guess  p,  see  if  it  works, 
and  then  do  some  revision. 

From  the  values  of  b  and  f  (for  all  i)  used  in  this 
paper,  we  can  compute  the  personnel  requirements,  shown  below. 


Manpower  Class 

1  2  3  4  5 


n 

286 

2156 

— 

713 

981 

537 

Rfl 

706 

20  34 

794 

2049 

1253 

A 

ut 

1264 

2540 

1028 

4093 

1790 

G 

1046 

1988 

1019 

2376 

991 

E 

Kfl 

722 

1316 

657 

1328 

540 

Kfl 

414 

655 

283 

- 

338 

249 

Table  12-1.  The  personnel  inventories,  pik 
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13.  Stage  Substitution 


The  astute  reader  will  notice  that  we  have  not  allowed 
substitution  between  stages.  Thus  a  type  (i,k)  person  is  not 
allowed  to  fill  a  type  (i  +  1,  j)  billet  or  a  type  (i-1,  j) 
billet. 

Our  formuation  actually  allows  such  assignments,  however 
we  are  assuming  that  they  net  out  to  zero.  That  is  for  each 
(i,k)  person  filling  a  stage  i+1  billet,  there  is  an  (i+l,k) 
person  filling  a  stage  i  billet. 

One  way  to  study  the  possibility  of  stage  substitution 
is  to  vary  the  length  of  the  intervals  (s^_^,  that 

define  the  stages.  However,  the  billet  requirement  are  pre¬ 
sumably  set  with  some  experience  level  in  mind.  Therefore 
major  shifts  in  the  definitions  of  stages  would  require  adjust¬ 
ments  in  the  billet  requirements. 


14.  Cost 


The  cost  of  any  allocation  can  be  calculated  using 


I  l 


ik  ik*k  ' 


where  c^k  is  the  annual  cost  of  an  (i,k)  person  and  w^kyk  =  zik 
is  the  number  of  (i,k)  people.  The  cost  data  c^k  should 
include  salary,  benefits,  retirement;  training,  promotion,  and 
recruitment  costs.  Notice  we  can  define  ck  as  £ k  cik  ik , 
and  then  rewrite  the  cost  as  £k  ckyk-  We  interpret  ck  as 
the  average  career  cost  of  each  accession  into  manpower 
class  k. 
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15.  Review 


We  have  constructed  our  model.  The  remainder  of  the 
paper  is  devoted  to  ways  of  calculating  particular  assign¬ 
ments  and  contains  a  sample  calculation. 

We  summarize  here  the  notations  and  definitions  pre¬ 
sented  to  this  point. 


manpower  class 


indexed  by  k  =  1,2, ...,K 


stages  of  a  career 


jobs 


indexed  by  i  =  1,2,...,  I 
defined  by  length  of  service 

(s^,  •••  •  ® j) 

indexed  by  j  =  1,2, . . .  ,  J 


people 

accessions 


indexed  by  stage  and  class  (i,k) 

yk,  the  number  of  people  entering 
class  k  per  year 


stage  inventory 


billets  requirements 


wikyk  is  the  ntunber  of  (i»k) 
people;  y^  is  the  accession  rate 
and  w^  is  obtained  from  the 
survivor  curve  and  the  stage 


definition.  Note 


Zik  "  wikyk* 


indexed  by  stage  and  job  (i,j); 
b. .  is  the  billet  requirements 


indexed  by  (i,j,k);  f£  is  the 


fraction  of  (i,j)  billets  that 
should  be  filled  from  class  k. 


billet  sharing 


target  allocation 


actual  allocation 


class  sums 


job  sums 


people  sharing 


cost 


t£j.  The  number  of  (i,k)  people 
desired  in  billet  (i,j) 


a£j  the  number  of  (i,k)  people 
in  billet  (i,j) 

and  b^  are  respectively 
the  number  assigned  to  billet  (i,j) 
and  the  requirement  for  billet 
(i#  j) 

zik  ^ik  are  resPectively 

the  number  of  (i,k)  people, 
and  the  target  for  (i,k)  people. 

gkj'  the  fraction  of  (i,k)  people 
assigned  to  billet  (i,j); 

akj  =  gkj  2ik* 

cik  is  the  annual  cost  of  an 

(i,k)  person.  The  cost  of  a  particular 

allocation  is  ^ 
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16.  Choice  of  an  Allocation 

Up  to  this  point  we  have  set  out  a  model  structure. 

The  next  task  is  to  devise  one  or  two  more  procedures  for 
selecting  allocations.  There  are  two  general  approaches  to 
this  task. 

On  approach  is  to  use  an  ad  hoc  rule.  For  example, 

we  could  fix  accessions  at  projected  rates  for  each  class. 

That  would  fix  the  variables  and  thus  z^k  = 

Next  we  could  specify  the  people  sharing  rules  g^ . .  That 

k  3 

would  give  us  the  allocation  aj^  *  gkjzik* 

The  second  approach  is  to  use  some  type  of  optimization 
to  select  a  policy.  This  optimization  is  either  based  on  a 
trade-off  between  the  cost  of  any  allocation  and  some  measure 
of  its  quality,  or  simply  a  measure  of  the  quality  of  the 
allocation. 
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17.  Quality  of  an  Allocation 

In  order  to  choose  an  allocation  using  optimization 
we  need  some  measure  of  che  quality  of  that  allocation. 

In  our  model  we  have  already  decided  to  focus  on  the 
discrepancy  between  people  inventory  and  targets  (z^k  -  PiJc) 
and  the  discrepancy  between  billet  assignments  and  requirements 

(xi j  "  bij) * 

We  have  selected  the  simplest  form  of  quality  measure 
that  makes  sense.  It  is  a  penalty  function  that  is  zero  for 
a  perfect  allocation  and  positive  for  others.  Thus  it  measures 
the  lack  of  quality.  The  penalty  function  has  three  properties. 

(1)  The  penalty  can  be  written  as 


l  l  hik(*ik>  +  l  l  *tj  <xij» 


where 

(2)  h..(2.v)  and  £.  .  (x.  .)  are  nonnegative,  convex  and  equal 

ik  ik  i]  13 

to  zero  if  z^  ■  p^  and  x^  =  b^  respectively. 

(3)  Parameters  and  ^  measure  a  unit  upper  and  lower 

percentage  error  in  meeting  billet  requirements 


if 

x.  . 

-  (i +  ■ 

l. , (x. .)  =  1 

ID 

ID  ID 

or 

x.  . 

ID 

*  (1  -  • 

Similar  parameters 

*ik 

and 

5 are  used  for  discrepan 

cies  between  actual  inventories  of  people  and  require¬ 
ments  . 

25 


f 


zik  -  (1  "  *ik>pik  • 
zik  *  (1  '  4 ik)pik  ’ 

Item  (1)  above  says  that  the  quality  measure  is 
separable,  item  (2)  that  is  nonnegative,  convex,  and  zero  if 
the  assignment  is  exactly  on  target.  The  third  item  requires 
more  discussion.  Item  (3)  is  designed  to  answer  the  question: 
How  do  we  compare  a  4  percent  shortfall  in  meeting  a  critical 
target  with  a  10  percent  short  fall  in  meeting  a  less  critical 
target?  Our  answer  to  this  question  is  to  take  a  single 
target  as  a  benchmark  and  to  say  arbitrarily  for  that  target 
that  a  certain  percent  over  and  a  certain  percent  under  the 
target  yields  an  error  of  one.  Then  for  any  other  target 
we  can  compare  the  seriousness  of  deviations  with  our  benchmark. 
We  say  a  fraction  0„  over  target  or  i|k  ^  under  target  is 
as  serious  as  the  deviations  we  have  established  for  our 
benchmark . 

For  example,  we  could  take  stage  3  pilots  as  our  bench¬ 
mark  (i  =  3,  j  =  5) ,  and  take  the  percentage  under  as  4% 

(i^35  =  0.04)  and  percentage  over  as  10%  (035  =  0.10)  as  the 
unit  serious  deviation  for  our  benchmark.  Now  take  any  other 
category,  for  example,  stage  4  general  warfare  billets 
(i  =  4,  j  =  2).  Then  we  ask  how  much  under  target  would  the 
assignment  to  these  billets  have  to  be  in  order  to  be  as 


hik(2ik>  = 


if 


or 
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serious  as  a  4%  shortage  of  stage  3  pilots.  In  this  way  we 
can  try  to  make  the  essential  judgments  about  the  trade-offs 
between  different  categories. 

We  shall  now  give  two  examples  of  penalty  functions 
and  then  give  some  practical  examples  of  how  such  a  criterion 
could  be  used. 
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18.  Piecewise  Convex  Penalties 


One  of  the  simplest  ways  to  construct  the  penalty  measure 
for  deviations  from  target  is  to  use  a  piecewise  convex  function. 
To  simplify  notation  we  should  consider  the  case  of  a  discrepancy 
between  the  inventory  of  people  and  desired  inventories  and 
we  shall  drop  the  (i,k)  subscripts. 

We  need  a  piecewise  linear  convex  function  h(z)  that 
satisfies 


(i) 

h  (z)  =  0 

if 

x  = 

p 

(ii) 

h  ( z)  =  1 

if 

X  * 

(1  +  4>)p 

(iii) 

h  ( z)  =  1 

if 

X  * 

(1  -  6  )p 

Such  a  function  is  shown  in  Figure  2.  This  function 
can  be  represented  in  several  ways.  One  of  the  simplest  is  as 
the  maximum  of  two  linear  functions 


h  (z) 


max 


(EL£) 

P 


This  functional  form  is  flexible  and  easy  to  manipulate. 
We  shall  see  below,  that  it  can  be  used  to  obtain  an  allocation 
using  a  linear  program.  Two  examples  of  this  are  contained 
in  the  appendix. 
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Quadratic  Penalty  Functions 

Quadratic  penalty  functions  are  also  quite  easy  to  use. 
Unfortunately,  we  cannot  get  a  quadratic  function  to  satisfy 
all  of  our  conditions.  For  example  let  us  write 

£.(x)  =  a(x  -  b)^  +  8  (x  -  b)  +  y. 


We  want  to  have 


(i) 

l  (x)  =0 

if 

x  = 

b 

(ii) 

l  (x)  =1 

if 

x  = 

(1  +  O) b 

(iii) 

Mx)  =1 

if 

X  = 

(1  -  <Jj) b 

(iv)  1 (x)  >0 

where  item  (iv)  implies  {,  (x)  has  its  minimum  at  x  =  b. 

These  four  conditions  cannot  be  met  by  a  quadratic 
function  unless  0  =  4';  that  is  unless  there  is  a  symmetry 
between  being  under  and  over  target.  In  the  symmetric  case 
we  can  write  l (x)  as 


i  (x) 


when  0  =  4». 


When  the  penalties  are  not  symmetric,  we  must  relax  one  of 
our  four  conditions.  If  we  relax  condition  (iv)  then  JC  (x) 
becomes 
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...  1  /  x  -  b\2  1  (6  -  Hi)  ,x_-b. 

h(x)  =  04)  \  E  j  "  7  W  (  b  *  * 

This  function  has  its  minimum  midway  between  (1  +  6)b  and 
(1  -  ip)b  and  it  has  a  value  of 

(e  -  »)2 
7&P 

at  the  minimum. 

This  approximation  is  useful  if  6  and  ip  are  similar. 
However,  the  approximation  becomes  much  worse  when  6  and  4 
are  much  different. 

Two  quadratic  penalty  models  are  described  in  the 
appendix. 
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20.  Cost  Quality  Trade-Off 

We  saw  above  that  any  particular  allocation  has  an 
annual  cost  and  also  a  penalty  associated  with  its  deviation 
from  targets  for  billet  requirements  and  personnel  inventories. 
By  placing  different  weights  on  the  cost  and  penalty  we  can 
obtain  a  family  of  objectives  that  will  lead  to  efficient 
allocations . 

The  cost-penalty  frontier  is  shown  below.  Notice 
the  minimum  cost  solution  is  zero,  since  we  would  not  have 
any  accessions;  i.e.  y  =  0,  thus  no  cost. 

We  write  our  objective  as 


As  the  parameter  X  increases  from  0  to  1  more  emphasis  is 
placed  on  the  cost  and  less  on  the  penalty.  At  X  =  0,  we 
minimize  the  penalty  and  ignore  the  cost;  at  X  =  1  we 
minimize  cost  and  ignore  the  penalty.  For  intermediate  values 
of  X  we  establish  a  trade-off  between  cost  and  penalty. 

This  curve  should  be  used  with  some  caution  since  the  costs 
are  expressed  in  real  units  (dollars  per  year) ,  but  the 
penalty  is  expressed  in  ad-hoc  units. 
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penalty 


minimum 

penalty 


APPENDIX  A 

LINEAR  PROGRAMMING  MODELS 


This  appendix  describes  two  linear  programming  approaches 
for  selecting  an  allocation.  The  discussion  is  brief  and  is 
intended  for  those  familiar  with  linear  programming.  The  first 
model  is  called  LPB  for  billet  sharing  and  the  second  is  called 
LPP  for  people  sharing. 

In  both  models  the  idea  is  to  choose  the  number  of 
people  assigned  to  billet  (i,1),  and  the  numbe  -  of  accessions 
in  manpower  class  k;  these  are  respectively  denoted  and 

y^.  The  allocation  a^  is  then  fixed  using  an  ad  hoc  rule. 

We  have  selected  two  rules  and  thus  get  two  different  linear 
programming  models. 

The  first  rule  is  for  sharing  billets 
Rule  B  (billet  share) 


(A- 1 ) 


a 


M 


f  1  v 

kjXij 


Recall  that  fj^  is  the  fraction  of  type  (i,j)  billets  that 
we  desire  to  have  satisfied  from  manpower  class  k.  Rule  B 
allocates  the  error  in  meeting  billet  requirements  among  man¬ 
power  classes  so  that  there  is  a  constant  percentage  error  in 
each  class's  allocation.  If  there  is  a  percentage  error  of 
n . ^  in  meeting  the  billet  (i,j)  requirement,  i.e. 


14 


n 


ij 


# 


then  the  billet  sharing  rule  means  that  the  percentage  error 
in  the  assignment  of  type  k  manpower  to  billet  (i,j)  will 
also  be 

The  people  sharing  rule  is 


Rule  P 


(A—  2 ) 


akj  ^kjzik 


where 


:ik  =  is  the  number  of  type  (i,k) 
type  (i,k)  people  are  allocated  among  different 
people  sharing  parameters  gjV.  If  there  is  an 
in  the  inventory  of  type  (i,k) 


people.  The 
jobs  using  the 
error  of  nik 


nik  = 


(zik  ~  Pjk> 

Pik 


then  the  percentage  error  in  the  assignment  of  type  (i,k) 
people  to  job  j  will  be  n^k  (recall  t£j  =  g^jP^k) •  The 
billet  sharing  rule  leads  to  a  larger  linear  program  than  the 
people  sharing  rule. 

The  constraints  of  the  linear  program  for  billet 
sharing  are  given  below.  The  objective  is  the  same  for  both 
programs;  it  will  be  discussed  after  we  examine  the  constraints. 


I 
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The  first  set  of  constraints  is  simply  the  conservation 
of  people.  There  will  be  w^kyk  type  (i,k)  people  and  under 
the  billet  sharing  rule  f^.x. .  of  them  will  be  assigned  to 

K]  1  j 

billet  (i,j). 

The  second  set  of  constraints  measures  the  difference 
between  actual  inventories  of  type  (i,k)  people  w^kyk  311(1 
the  requirement  p^.  The  objective  is  selected  so  that  the 
solution  to  the  linear  program  will  have  eik  positive  and 
dik  equal  to  zero  if  w^kyk  exceed  pik  and  e^  equal  to 


zero  and 

dik 

positive  if 

wik*k  is 

less 

than  pik . 

The 

e . .  and 
ik 

dik 

thus  measure 

the  excess 

and 

deficit  in 

the 

type  (i,k)  people  account. 

The  variables  m. .  and  n. .  play  a  similar  excess 
and  deficit  role  in  the  type  (i,j)  billet  account. 

The  first  set  of  constraints  for  the  LP-P  program  is 


(A-3) 


xij 


K 

-  I 


kix  gkjWik*k 


=  0 


for  all  (i,j)  . 


This  again  is  a  conservation  constraint.  There  are 
wik^k  (i'k)  people  and  gj^  of  them  are  assigned  to 

job  j.  Thus  these  constraints  assure  that  the  assignments 
to  billets,  j »  are  consistent  with  actual  manpower  available. 

The  second  and  third  sets  of  constraints  in  LP-P 
are  identical  with  the  second  and  third  sets  of  constraints 
in  LP-B.  They  define  the  excess  and  deficit  in  people  and 
billet  accounts. 
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The  form  of  the  conservation  constraint  in  LP-P  allow 
for  some  simplifications.  Notice  that  if  all  the  yk  are 
nonnegative  then  (A-3)  implies  that  all  the  x^j  will  be 
nonnegative.  Thus  we  can  use  (A- 3)  to  eliminate  from 

the  problem.  With  this  simplification  the  constraints  of 
LP-P  Are  shown  in  Table  A-2. 

The  objective  for  both  linear  programs  will  typically 
be  a  compromise  between  minimizing  costs  and  minimizing  the 
penalty  that  measures  our  departure  from  people  and  billet 
targets.  The  cost  objective  is 


(A-4) 


K 

l 

k=l 


c,  y, 
kJk 


The  penalty  measure  is 
I  J 


(A- 5) 

l 

i=l  j 

h  wunij +  "ij-ij 

where 

(i) 

"ij  -  ,eijbiirl 

(ii) 

(A-6) 

(iii) 

Yik  ■  ^tkPik1"1 

(iv) 

°ik  *  (6ikpik) 
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Recall  that  6,  <f>,  and  <S  are  the  percentage  errors  that 

lead  to  a  unit  loss. 

A  balance  between  these  two  objectives  is  obtained  by 
talcing  a  weighted  combination. 


K 

(A-7)  A  l  c.y.  +  (1  -  A) 
k=l  K  K 


x 


I  J 

{  V  y  IT .  ,n,  .  +  y.  ,m.  , 

iii  j-i 


I  K 

ill  kh  ''ikeik  +  "ikdik’ 


The  problem  LP-B  has  I  *  (J  +  2R)  constraints  and 
K  +  I  x  (2K  +  3J)  variables.  For  our  example  with  1=6, 

J  =  7  and  K  =  5  this  means  102  constraints  and  191  variables. 

The  problem  LP-P  has  I  x  (J  +  K)  constraints  and  K  +  1  x  (2K +  2j) 
variables.  For  our  example  this  works  out  to  72  constraints 
and  149  variables.  LP-P  is  particularly  easy  to  solve  since 
one  can  always  get  a  reasonable  first  basic  solution  by  guessinq 
the  y^  and  then  choosing  the  eik'  mij'  and  nij  to 

satisfy  the  constraints. 

The  duals  of  both  LP-B  and  LP-P  appear  to  be  easier  to 


solve.  The  dual  variables  for  each  problem  and  each  set  of 
constraints  is  shown  in  Tables  (A-l)  and  (A-2)  respectively. 
The  dual  program  for  LP-B  is 


(A- 8)  maximize 


I  K 


I  J 


iii  Ji  ViltPilt  +  Ji  jli  qi5bli 


subject  to 


^  uik£«  +  ■’ij  i 0 


for  all  (i,j) 


I  I 

■  uikwik +  vikwik  i  Xck  tor  a11  * 


-  <1“x)'Yik  £  vik  £  (1  -  x)oik  for  all  (i,k) 

-  (l-X)ir..  <  q..  <_  (1  -  X)  y .  .  for  all  (i,j) 

lj  lj  ij 

Dual  of  LP-B 


The  dual  of  LP-B  has  K  +  I  x  J  constraints,  I  *  (2K  +  J) 
variables,  and  I  *  {K  +  J)  of  these  variables  have  upper  and 
lower  bounds.  In  our  example  this  works  out  to  47  constraints, 
102  variables,  and  72  variables  subject  to  upper  and  lower  bound 
constraints. 

The  dual  of  LP-P  is  shown  below. 


(A-9) 


This  problem  has  K  constraints  and  I  *  (J  +  K)  variables  with 
upper  and  lower  bounds.  In  our  example,  this  would  be  5  con¬ 
straints  and  72  variables  with  upper  and  lower  bounds.  This 
problem  may  be  quite  easy  to  solve. 
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APPENDIX  B 


UNCONSTRAINED  QUADRATIC  MODELS 


This  appendix  outlines  two  methods  for  selecting  an  allo¬ 
cation  using  the  quadratic  penalty  functions.  The  two  models  are 
similar  to  the  two  linear  programming  models;  they  use  the  billet 
sharing  or  people  sharing  rules  to  go  from  an  aggregate  problem 
to  a  detailed  allocation. 

The  models  do  not  have  any  inequality  contraints.  This 
is  to  insure  that  the  solution  can  be  obtained  quickly.  A  full 
model  with  inequality  constraints  might  appear  to  be  more  appro¬ 
priate,  however,  we  must  recall  our  ultimate  objective  is  to  cal¬ 
culate  allocations  in  a  relatively  simple  way.  The  model  is  not 
built  on  exact  premises  and  it  does  not  use  precise  data.  It  is  there¬ 
fore  not  terribly  important  to  be  exact  in  choosing  an  allocation. 

We  hope  to  have  an  easy  and  consistent  way  of  choosing  allocations 
so  we  can  compare  the  effects  of  changing  assumptions  on  the 
allocations. 

Both  models  calculate  the  billet  assignments  x^  and 
the  accession  rules  y^.  The  actual  allocation  is  determined  by 
the  billet  sharing  rule  in  model  UQ-B 
Rule  B. 


(B-l) 


akj  =  fkjxij 
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In  model  UQ-P,  the  people  sharing  rule  is  used 
Rule  P. 


(B-2) 


*kj  ~  gkjwikyk 


These  rules  are  discussed  in  Appendix  A. 

The  objective  in  these  quadratic  models  is  based  on 
symmetric  penalties;  i.e.  9^  =  and  4^  =  6^^.  The 

objective  is  a  combination  of  two  terms:  the  cost  term  and 
a  penalty  term 

K 

( B—  3)  A  l  ckyk  +  (1  -  X) 


k=l 


I  J 


I  K 


(  &  E  ^ii^ii  '  bii)  +  E  E  Yik(wikYk  “  Pik}  } 

i=l  j=l  13  13  13  i=l  k=l  1K  xK  K  1K 


where 


(B-4) 


(i)  vi..  =  (0..b..) 
13  13  13 


(ii)  Yik  -  <*ikPik> 


-2 


-2 


In  the  UQ-B  model  we  minimize  ( B—  3 )  subject  to  the  conservation 
constants. 


(B-5) 


■  i 

fkjxij  “  wikyk  =  0  for  each  ( i , k ) 
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The  people  sharing  model,  UQ-P,  minimizes  (B-3)  subject 


to  the  constraints 


(B-6) 


x. 


K 

-  I 

k=l 


gkjwikyk 


=  o 


In  UQ-P  we  can  substitute  (B-6)  directly  into  (B-3)  and  thus 
solve  a  completely  unconstrained  problem.  In  UQ-B  we  minimize 
(B-3)  subject  to  (B-5) . 


r 

» 

APPENDIX  C 
SAMPLE  SOLUTION 

This  section  contains  some  sample  solutions  using  the 
unconstrained  minimization  described  in  Appendix  B. 

The  data,  b^j,  f£j,  wifc'  P^'  are  described  in  the  paper. 
The  p^k  and  g£,  are  calculated  in  the  manner  suggested  in 
Section  12. 

The  solutions  are  based  on  data  that  is  largely  subjective 
and  does  not,  in  any  way,  pretend  to  capture  the  situation  that 
exists  in  the  Navy.  The  intent  of  this  section  is  to  demonstrate 
the  feasibility  of  the  scheme. 

The  target  errors  are  expressed  in  percentage  terms.  The 
10,000%  target  error  indicates  that  we  do  not  care  very  much  about 
meeting  billit  or  personnel  targets  in  the  first  two  stages. 


Job 

Stage"'>sv^ 

1 

2 

3 

4 

5 

6 

7 

ENS 

10000 

10000 

10000 

10000 

10000 

10000 

10000 

LTJG 

10000 

10000 

1U000 

10000 

10000 

10000 

10000 

LT 

10 

8 

6 

4 

4 

5 

6 

LTCDR 

8 

6.4 

4.8 

3.2 

3.2 

4 

4.8 

CDR 

7 

5.6 

4.2 

2.8 

2.8 

3.5 

4.2 

CAPT 

6 

4.8 

3.6 

2.4 

2.4 

3 

3.6 

Table  C-l .  The  inputs  9^  in  %;  for  billet  requirements 


Job 

Stage _ 1 


2  3  4  3 


ENS 

0 

10000 

LTJG 

10000 

10000 

LT 

LTCDR 

10 

8 

7.5 

6 

CDR 

CAPT 

7 

5.25 

6 

4 . 5 

10000 

10000 

1000C 

10000 

10000 

10000 

4.5 

5 

6 

3.6 

4 

4.8 

3.15 

3.5 

4.2 

2.7 

3 

3.6 

ik 


Table  C-2/  The  inputs  $ 


in  % ;  for  manpower  inventories 
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The  billet  assignments  and  the  actual  percentage  errors 


using  the  people  share  rule  are 


STAGE/JOB 

1 

2 

3 

4 

5 

6 

7 

ENS 

1444 

0 

2231 

1000 

1791 

972 

0 

LTJG 

1498 

0 

1989 

895 

1669 

908 

1 

LT 

2164 

383 

2380 

19  76 

2021 

992 

399 

LTCDR 

1242 

444 

1246 

586 

977 

395 

452 

CDR 

1108 

482 

1000 

489 

5  78 

63 

884 

CAPT 

540 

529 

308 

150 

0 

0 

510 

Table  C-3.  The  billet  assignment, 


STAGE/JOB 

1 

2 

3 

4 

5 

6 

7 

ENS 

"253 

0 

6 

“45 

"86 

"87 

0 

LTJG 

"  47 

0 

6 

"22 

"16 

25 

20 

LT 

20 

1 

"14 

“27 

43 

29 

37 

LTCDR 

17 

26 

18 

28 

45 

28 

38 

CDR 

7 

3 

8 

2 

15 

"15 

2 

CAPT 

9 

8 

4 

0 

0 

0 

"32 

Table  04.  The  percentage  error  in  meeting  billet  requirements 


The  accessions  to  the  five  classes  are: 

Class  1 _ 2 _ 3 _ 4 _ 5 

|  645  1162  524  926  508 

Table  C-5.  Annual  rate  of  accessions,  yk 
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These  accessions  produce  the  actual  personnel  inventories 
and  the  precentage  error  in  meeting  inventories  are  given  below. 


Class 

1 

2 

3 

4 

5 

ENS 

1274 

2296 

1035 

1829 

1003 

LTJG 

1192 

2149 

969 

1712 

938 

LT 

1613 

2906 

1311 

2316 

1269 

LTCDR 

909 

1637 

738 

1304 

715 

CDR 

789 

1421 

641 

1132 

621 

CAPT 

349 

629 

284 

501 

275 

Table  C-6.  Personnel  inventory  by  class  and  stage,  • 


Table  C-7.  Percentage  errors  in  meeting  personnel  inventory  budgets. 


The  solution  using  the  billet  sharing  rule  produced  similar 
results.  The  billet  assignments  and  %  error  in  billet  assignments 
were 
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T 


STAGE/JOB  123456  7 


ENS 

1801 

0 

2074 

940 

1388 

1015 

0 

LTJG 

1685 

0 

1941 

879 

1298 

949 

1 

LT 

2280 

405 

2423 

1109 

1399 

1021 

499 

LTCDR 

1284 

370 

1294 

596 

6  76 

493 

432 

CDR 

1115 

496 

1036 

483 

101 

64 

1173 

CAPT 

493 

554 

291 

147 

0 

0 

492 

Table  C-8 


STAGE/ JOB 

1 

2 

3 

4 

5 

6 

7 

ENS 

"340 

0 

1 

"37 

"44 

"95 

0 

LTJG 

67 

0 

3 

“20 

35 

22 

0 

LT 

26 

7 

"17 

"31 

61 

27 

21 

LTCDR 

14 

38 

12 

26 

62 

10 

40 

CDR 

8 

5 

"12 

4 

85 

"17 

"30 

CAPT 

17 

"13 

9 

2 

0 

9 

"27 

Table  C-9 


The  accessions  using  billet  sharing  were 


Class 


2  3  4 

1187  530  748 

Table  C-10 


5 

550 
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These  accession  rates  give  rise  to  the  personnel 


inventory  and  percentage  error  in  personnel  levels 


Class 

1 

2 

3 

4 

5 

ENS 

1261 

2345 

1048 

1478 

1087 

LTJG 

2194 

981 

1383 

1017 

LT 

1596 

2968 

1326 

1871 

1376 

LTCDR 

899 

1672 

747 

1054 

775 

CDR 

1451 

649 

915 

673 

CAPT 

|  345 

642 

287 

405 

29  8 

Table 

C-ll 

Class 

1 

2 

3 

4 

5 

ENS 

"340 

9 

47 

51 

~102 

LTJG 

67 

8 

23 

33 

19 

LT 

26 

“17 

29 

54 

23 

LTCDR 

14 

10 

1 

31 

24 

CDR 

17 

2 

2 

20 

"20 

Table  C-12 


The  analysis  of  this  output  should  be  directed  toward 
constructive  changes  in  the  input  data.  Can  we  shift  some 
billet  assignments  from  LT  to  LTJG  and  LTCDR?  Should  we 
tighten  up  more  on  the  pilot  inventory?  Should  we  change 
the  sharing  rules?  Is  it  possible  to  alter  the  survivor 
fractions,  and  thereby  improve  the  solution? 


We  see  that  there  are  a  number  of  potential  questions 


we  can  answer.  The  model  is  flexible  and  the  calculations 
are  rapid.  It  should  be  an  excellent  tool  for  analyzing 
manpowe r  policy. 
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FILE  STRUCTURE 


Component 

Number 

Dimension 

Symbol 

Description 

1 

(K,I,J) 

G 

People  sharing  array.  +/G  is 

an  ( I , K)  element  matrix  of  ones. 

G  is  nonnegative:  gj^  is  the 
fraction  of  type  (i,k)  people 
assigned  to  job(i,j) 

2 

(I,J) 

B 

The  billet  requirements,  b^ 
is  the  number  of  type(i,j)  jobs 

to  be  filled. 

3 

(I,J) 

0 

The  range  for  a  unit  error  in 

overfilling  billets.  If 

x.  .  =  (1  +  0.01  .  .)b.  .  there 
lj  lj  xj 

is  a  unit  error.  0^.  >0. 
Expressed  as  %. 

4 

(I,K) 

The  range  for  a  unit  error  in 

overfilling  personnel  inventories 

2ik  =  <ll0-0Uik)pik  a  unit 

error  is  counted.  >}>..  >0. 

lk 

Expressed  as  %. 

5 

(K,  I , J) 

F 

The  job  sharing  array.  F  is 

nonnegative  and  +/[1]F  is 

an  (I,J)  matrix  of  ones:  f t  . 

kj 

is  the  fraction  of  type(i,j) 
jobs  to  be  filled  by  type(i,k) 
people. 

6 

(If  K) 

P 

The  desired  manpower  inventory: 
p^k  is  the  number  type(i,k) 
people  desired. 
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FILE  STRUCTURE  CONT 


Component 

Number 


Dimensional  Symbol 


(N,K) 


(I,K) 


(I,K) 


(K ,  I ,  J) 


(I,J) 


(K  ,  I ,  J) 


Description 


The  survivor  fractions  for 
each  of  the  K  classes .  N 
is  the  maximum  LOS. 


The  last  year  in  each  stage, 
Increasing  positive  numbers 
with  Sj  =  N . 


The  element  w., 
lk 


is  the  number 


of  years  a  class  k  input  expects 
to  spend  in  stage  i. 


The  calculated  accessions  rate; 
y^  is  the  number  of  accessions 
in  manpower  class  k. 


The  calculated  inventory  of 
each  type  of  person.  There  are 
zik  type  people. 


The  calculated  allocation:  a.1 
is  the  number  of  type  (i,k) 
people  filling  type(i,j)  jobs. 


The  calculated  billets  filled. 
There  are  x^  people  in  job(i,j). 


The  target  allocation.  We  want 
type(i,k)  people  in  type(i,j) 
jobs.  If  things  are  consistent 
+/T  equals  P  and  +/[1]T  equals  B. 


The  range  for  a  unit  error  in 

under  filling  billets.  If 

x^  =  (l-i|'^)b^j  then  an  error  of  1 

is  counted .  u>.  .  >  0 . 

ID  - 


FILE  STRUCTRE  CONT 


Component 

Number 

Dimensional 

Symbol 

Description 

16 

(I,K) 

6 

The  range  for  a  unit  error  in 
underfilling  personnel  inven- 
tories .  If  *ik  -  <l-«ik>Pik 

then  a  unit  error  is  counted. 

6  ,k  >  °. 

17 

(K, I , J) 

c 

The  element  c£_.  is  the  total 
cost  of  having  a  type  (i,k) 
person  in  job(i,j)  for  one 

year . 

J 


FUNCTIONS 


Function 

Uses 

Computes 

Syntax  and  Description 

BILQUAD 

B,  6,  <f»,  F, 

P,  W 

Y,  X,  Z,  A 

BILQUAD  ’FILENAME* 

Calculates  "optimal"  acces¬ 
sions  y,  personnel  inventory 
(2),  and  allocation  (A), 
given  the  data,  objective, 
and  the  billet  share  rule  (F) . 

BILSIM 

F,Y,X 

Z ,  A 

BILSIM* FI LEN AMS’ 

Given  Y  and  X  (which  pre¬ 
sumably  satisfy  the  conser¬ 
vation  constraint) ,  calculates 
allocation  (A)  and  personnel 
inventory  (Z),  using  the 
billet  sharinq  rule. 

PGREC 

P.  G 

B,  F,  T 

—  - 

PGREC  ' FILENAMS ' 

reconciles  the  file  elements 

B,  G,  and  T  with  P  and  G 

BFREC 

B,  F 

■ 

BFREC’ FILENAMS' 

reconciles  the  file  elements 

P,  G,  and  T  with  the  file 

elements  B  and  F  in  the 

same  file. 

PEOQUAD 

G,  B,  P,  0, 

F,  P,  W 

PEOQUAD  ’FILENAMS’ 

calculates  optimal  acces¬ 
sions  (y) ,  personnel  in¬ 
ventory  (z),  billet  staffing 
(x) ,  and  allocation  (A) , 
given  the  data,  and  use  of 
the  people  share  rule. 
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FUNCTIONS  CONT. 


FUNCTION 

USES  ~~ 
COMPONENTS 

CALCULATES 

COMPONENT 

SYNTAX  AND  DESCRIPTION 

PEOSIM 

G,  Y,  W 

Z,  X,  A 

PEOSIM 'FILENAMS' 
calculates  inventory  z, 
staffing  x  and  allocation 

A,  where  accessions  and 

people  share  rule  are  given. 

WAITS 

a,  S 

W 

WAITS  'FILENAMS' 

Given  the  survivor  frac¬ 
tions  (a)  and  stage  defini¬ 
tions  (s) ,  calculates  w 
the  expected  waiting  time 

in  each  state. 

PEOPRONT 

P,  z 

lOOx(D-z) *P 

PEOPRONT  'FILENAMS' 

calculates  the  percentage 
error  in  inventory.  Compare 
with  ip . 

BILPRONT 

B,  X 

100x(b-x)tB 

1 

BILPRONT  'FILENAMS' 

calculates  the  percentage 
error  in  meeting  billet 

requirements.  Compare 

with  0 . 

ROUND 

ROUND  XXX 

takes  any  array  and  rounds 
elements  to  integers. 

REFERENCES 


[1]  Grinold,  R.  C. ,  and  K.  T.  Marshall.  Manpower  Planninq, 

North-Holland  (1977) . 

[2]  Holt,  C.C.,  F.  Modigliani,  J.  F.  Muth ,  and  H.A.  Simon, 

Planning  Production  Inventories  and  Work  Force, 
Prentice  Hall  (19(50) . 


INITIAL  DISTRIBUTION  LIST 

No.  of  Copie 


Defense  Documentation  Center  2 

Cameron  Station 
Alexandria,  VA  22314 

Library,  Code  0142  2 

Naval  Postgraduate  School 
Monterey,  Ca .  93940 

Library,  Code  55  1 

Naval  Postgraduate  School 
Monterey,  Ca.  93940 

Dean  of  Research  1 

Code  012 

Naval  Postgraduate  School 
Monterey,  CA  93940 

Director,  1 


Total  Force  Management  Control/Analysis  Div. 

(OP-10) 

Office  of  DCNO(MPT) 

Department  of  the  Navy 
Washington,  D.C.  20370 

Director, 

Total  Force  Planning  Div.  (OP-11)  1 

Office  of  DCNO(MPT) 

Department  of  the  Navy 
Washington,  D.C.  20370 

Director,  1 

Total  Force  Program  Div.  (OP -12) 

Office  of  DCNO(MPT) 

Department  of  the  Navy 
Washington,  D.C.  20370 

Director,  1 

Military  Personnel/Training  Division  (OP-13) 

Department  of  the  Navy,  Office  of  DCNO(MPT) 

Washington,  D.C.  20370 

Dr.  A.  Charnes  1 

Dept,  of  Math. 

University  of  Texas 
Austin,  TX  78712 


59 


r 


No.  O 


Prof.  M.  Cooper 

Dept,  of  Ind.  Eng.  and  CR 

Inst,  of  Technology 

Southern  Methodist  University 

Dallas,  Texas  75275 

Mrs.  Mary  Snavely  Dixon 
Office  of  DASN  (Manpower) 
Department  of  the  Navy 
Washington,  D.C.  20350 

LCDR  Lee  Gunn  (OP-132C1) 

Surface  Officer  Program  Manager 
Office  of  DCNO(MPT) 

Arlington  Annex 
Washington,  D.C.  20370 

LCDR  Gary  Johnson  (OP-136C1) 
Office  of  DCNO(MPT) 

Arlington  Annex 
Washington,  D.C.  20370 

Dr.  Robert  F.  Lockman 
Center  for  Naval  Analyses 
1401  Wilson  Blvd 
Arlington,  VA  22209 

Mr.  N.  K.  Malehorn  (OP-102) 

Office  of  DCNO(MPT) 

Arlington  Annex 
Washington,  D.C.  20370 

Professor  James  G.  March 
School  of  Business  Administration 
Stanford  University 
Stanford,  CA  94305 

Dr.  Kneale  T.  Marshall  (OP-01T) 
Office  of  DCNO(MPT) 

Arlington  Annex 
Washington,  D.C.  20370 

CAPT  Rodney  B.  McDaniel  (OP-110) 
Office  of  DCNO(MPT) 

Arlington  Annex 
Washington,  D.C.  20370 


R.  J.  Niehaus 

3300  N.  Agindon  St., 

Arlington,  VA  22207 


60 


l 


No.  of  Copies 

LCDR  Daniel  Parker  (OP-130D2)  1 

Office  of  DCNO(MPT) 

Arlington  Annex 
Washington,  D.C.  20370 

Alfred  S.  Rhode,  Vice  President  1 

Information  Spectrum  Inc., 

1745  South  Jefferson  Davis  Highway 
Suite  401 

Arlington,  VA  22202 

Dr.  Bernard  D.  Rostker  1 

PDASN  (M&RA) 

Department  of  the  Navy 
Washington,  D.C.  20350 

Joe  Silverman  (Code  303)  1 

Navy  Personnel  RsD  Center 
San  Diego,  CA  92152 

Prof.  Robert  E.  Stanford  1 

Dept,  of  Management 
Auburn  University 
Auburn,  AL  36830 

Prof.  R.  Grinold  10 

Dept,  of  Business  Administration 
University  of  California 
Berkeley,  CA  94720 

Dr.  T.  C.  Varley  (Code  434)  1 

Office  of  Naval  Research 
800  N.  Quincy  St. 

Arlington,  VA  22217 

Naval  Postgraduate  School 
Monterey,  Ca .  93940 


Attn:  J.  Arima,  Code  54  1 

C.  Jones,  Code  54  1 

P.  Carrick,  Code  54  1 

R.  Elster,  Code  54  1 

G.  Thomas,  Code  54  1 

L.  Cavaiola,  Code  55  1 

D.  Gaver,  Code  55  1 

P.  Jacobs,  Code  55  1 

P.  Lewis,  Code  55  1 

P.  Milch,  Code  55  20 

R.  Stampfel,  Code  55  1 

M.  Sovereign,  Code  55  1 


61 


